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Abstract 

We provide simple examples of closed-form Gaussian wavepacket solutions of the free-particle 
Schrodinger equation in one dimension which exhibit the most general form of the time- 
dependent spread in position, namely (Axt) 2 = (Axo) 2 + At + (Apo) 2 t 2 /m 2 , where A = 
{(x — (x)o)(p — {p)o) + (p — {p)o)(x — (£)o))o contains information on the position-momentum cor- 
relation structure of the initial wave packet. We exhibit straightforward examples corresponding 
to squeezed states, as well as quasi-classical cases, for which A < so that the position spread can 
(at least initially) decrease in time because of such correlations. We discuss how the initial corre- 
lations in these examples can be dynamically generated (at least conceptually) in various bound 
state systems. Finally, we focus on providing different ways of visualizing the x — p correlations 
present in these cases, including the time-dependent distribution of kinetic energy and the use of 
the Wigner quasi-probability distribution. We discuss similar results, both for the time-dependent 
Axt and special correlated solutions, for the case of a particle subject to a uniform force. 
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I. INTRODUCTION 



The study of time-dependent solutions of the one-dimensional Schrodinger equation is 
a frequent topic in many undergraduate textbooks on quantum mechanics. The problem 
of a Gaussian or minimum-uncertainty wavepacket solution for the case of a free particle 
(defined more specifically below) is the most typical example cited, often being worked out 
in detail, or at least explored in problems The emphasis is often on the time- dependent 
position spread for such solutions, typically written in the forms 

(Ax,) 2 = (Ax ) 2 + (£j 2 \ = ( Ax f + ^!!! (1) 

where the spreading time or coherence time can be defined by to = mAxo/ Ap^. Textbooks 
rightly point out the essentially classical nature of much of this result, explained by the fact 
that the higher momentum components of the wave packet outpace the slower ones, giving 
a position-spread which eventually increases linearly with time as Axt ~ Avot, where Avq 
is identified with Apo/m. 

The form of the expression for Axt in Eqn. (0) is a special case of the most general possible 
form of the time-dependent spatial width of a one-dimensional wave packet solution of the 
free-particle Schrodinger equation which is well-known in the pedagogical literature 
but seemingly found in many fewer textbooks |l0|. This general case can be written in the 
form 

(A^) 2 = (Ax ) 2 + ((x - (x) )(p - (p> ) + (p- (p) )(x - (x) )) - + { -^f- ( 2 ) 

where the coefficient of the term linear in t measures a non-trivial correlation between the 
momentum- and position-dependence of the initial wave packet. While such correlations 
are initially not present in the standard Gaussian wave packet example routinely used in 
textbook analyses, which therefore gives rise to the simpler form in Eqn. (0), a non- vanishing 
x—p correlation does develop for later times as has been discussed in at least one well-known 
text and several pedagogical articles 

For wave packets which are constructed in such a way that large momentum components 
(p > (p)o) ar e initially preferentially located in the 'back' of the packet (x < (x)q), the initial 
correlation can, in fact, be negative leading to time-dependent wave packets which initially 
shrink in size, while the long-time behavior of any ID free particle wave packet is indeed 



always dominated by the quadratic term in Eqn. (J2J), consistent with standard semi-classical 
arguments. (We stress that we will consider here only localized wave packets which are 
square-integrable, for which the evaluation of Ax t and Ap t is possible, and not pure plane 
wave states nor the special non-spreading, free-particle solutions discovered by Berry and 
Balazs Q.) 

For the standard Gaussian or minimum uncertainty wave packet used in most textbook 
examples, and in fact for any initial wave packet of the form ip(x, 0) = R(x) exp(ip Q (x—xo) / K) 
where R(x) is a real function, this initial correlation vanishes and the more familiar special 
case of Axt in Eqn. ((T} results, leading many students to believe that it is the most general 
result possible. It is, however, very straightforward to construct initial quantum states con- 
sisting of simple Gaussian wave functions, such as squeezed states or linear combination of 
Gaussians, which have the required initial position- momentum correlations 'built in', and 
which therefore exhibit the general form in Eqn. (J2J), including examples where the position- 
space wave packet can initially shrink in width. Since these examples can be analyzed with 
little or no more mathematical difficulty than the standard minimum-uncertainty cases com- 
monly considered in textbooks [ifl, we will focus on providing two such examples below. We 
will, however, also emphasize the utility of different ways of visualizing the time-dependent 
position- momentum correlations suggested by the form in Eqn. (0). 

The derivation of Eqn. (J2J) has been most often discussed Q|>la| using the evaluation of 
the time-dependence of expectation values described by 

using the free particle Hamiltonian, H = p 2 /2m, or related matrix methods j^J; since we 
are interested only in expectation values of operators (A = x or p) which are themselves 
independent of time, there is no additional (dA/dt) term in Eqn. In the next section, 
we derive the necessary time-dependent expectation values of powers of position and mo- 
mentum in a complementary way, using very general momentum-space ideas. (Identical 
methods can then also be used to evaluate the general form of Ax t for the related case of 
uniform acceleration, which we discuss in Appendix El) Then in Sec. IIHI we briefly review 
the special case of the minimum-uncertainty Gaussian wave packet (to establish notation) 
focusing on the introduction of useful tools to help visualize possible correlations between 
position and momentum in free particle wave packets, especially the direct visualization of 
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the real/imaginary parts of ip(x,t), the time-dependent spatial distribution of kinetic en- 
ergy, as well as the Wigner quasi-probability distribution. Then, in Sec. IIVI we exhibit two 
cases of correlated wave packets with the general form of Ax t in Eqn. (J2J), which are sim- 
ple extensions of these standard results. A similar example, involving squeezed states, has 
been discussed in Ref. |l4|, but we will focus here on understanding the detailed position- 
momentum correlations which give rise to the term linear in t in Eqn. (J2J), especially using 
the techniques outlined in Sec. IHII for their visualization. Finally, we make some concluding 
remarks as well as noting in an Appendix that very similar results (both for the general 
form of the time-dependent Axt and for the exemplary cases studied here) can be obtained 
for the Schrodinger equation corresponding to the case of constant acceleration. 



II. TIME-DEPENDENT Ax t USING MOMENTUM-SPACE WAVEFUNCTIONS 

While the general result for the free-particle Ax t is most often obtained using formal 
methods involving the time-dependence of expectation values as in Eqn. (jHJ), one can also 
evaluate time-dependent powers of position and momentum for a free particle in terms of 
the initial wave packet quite generally in terms of the momentum-space description of the 
quantum state, namely <p(p,t), obtaining the same results, in a manner which is nicely 
complementary to more standard analyses. Depending on the ordering of topics in a given 
quantum mechanics course syllabus, this discussion might well be applicable and under- 
standable earlier in the curriculum than the more formal method. 

In this approach, the most general momentum-space wave function which solves the free- 
particle time- dependent Schrodinger equation 

t) = H<f>{p, t) = E(f>{p, t) = ih-<P(p, t) , (4) 

can be written in the form 

<P(p,t) = Mp)e- ip2t/2mh (5) 

with <p(p, 0) = <j>o(p) being the initial state wavefunction. The t-dependent expectation 
values for powers of momentum are trivial since 

/+oo 
p\Mp)\ 2d P= (P)o (6) 
-oo 

/+oo 
p 2 \Mp)\ 2 d P =(p 2 )o (7) 
-oo 
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so that 

(A^) 2 = (p 2 ) t - (p) 2 = (p\ - (p)l = (A Po ) 2 (8) 

as expected for a free-particle solution for which \<f>(p, t)\ 2 = \cf> (p)\ 2 is independent of time. 

In this representation, the position operator is given by the non-trivial form x = %h{dj dp), 
and the time- dependent expectation value of position can be written as 



-+oo 

(x) t = I {<P(p,t)]*x{<P(p,t)]dp 

-oo 

"+oo / o \ 

<P*Jp) e +ip2t/2mh iih^-\ <j) (p) e- ip2t/2mh 
oo L J \ op J L 

+00 / o \ . n+00 



dp 



[0o (P)] ( «^ J [0o (p)] dp+ — I p\Mp)\ 2 dp 

= (^)o + £(p)o (9) 

which is consistent with Ehrenfest's theorem for the essentially classical behavior of (x) t . 
The same formalism can be used to evaluate (x 2 )t and gives 

t t 2 
(x 2 ) t = (x 2 ) + —(xp + px)o + (p 2 )o^ (10) 
m m z 

where one can use the general representation-independent commutation relation [x,p] = ih 
to simplify the answer to this form. The symmetric combination of position and momentum 
operators, written here as (xp + px) , which is obviously Hermitian, guarantees that this 
expression is manifestly real. (Discussions in textbooks on symmetrizing products of non- 
commuting operators abound, but such results are seldom put into the context of being 
useful or natural in specific calculations, as is apparent in their use here.) 

Combining Eqns. Q and (fTU|) then gives the most general form for the time- dependent 
spread in position to be 



(A*,) 2 = (x 2 ) t - (x)\ 



t 2 



(x ) H (xp + px) + (p ) — - (x) H (p) 

m m z J \ m 

= (Ax ) 2 + ((xp + px) -2(x) (p) )-+ { -^^ 

= (Ax ) 2 + ((x - (x) )(p - (p) ) +(p- (p) )(x - (x) )) - + . (11) 

We have rewritten the term linear in t in a form which stresses that it is a correlation 
between x and p, similar in form to related classical quantities such as the covariance in 



probability and statistics. Recall that for two classical quantities, A and B, described by a 
joint probability distribution, the covariance is defined as 



cov(A, B) = ((A - (A)) (B - (B))) = (AB) - (A)(B) . (12) 

As we will see in the next section, there is no initial correlation for the familiar minimum- 
uncertainty Gaussian wave packets. However, for simple variations on the standard example, 
as in Sec. IIV[ we will find non-vanishing correlations, which we can visualize with the 
methods in Sec. Mil 

We stress that the notion of a time-dependent correlation between x and p at arbitrary 



times (t > 0) can be easily genera 
covariance for these two variables 



ized from these results, and we can define a generalized 
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12j (or any two operators, A, B) as 



cov{x,p;t) = -((x- (x) t )(p- (p)t) + {p- (P)t){x- (x) t )) t (13) 



where the additional factor of 1/2 accounts for the necessarily symmetric combination which 
appears, compared to the classical definition. One can then speak of a time- dependent 
correlation coefficient defined by 

P{x,P]t) = — -— 14 

Ax t ■ Apt 



in analogy with related quantities from statistics. This correlation can be shown 
satisfy the inequality 



121 to 



' h - 2 



which vanishes for the standard minimum- uncertainty Gaussian at t = 0, but which is 
non-zero for later times, as we will see below. 



III. STANDARD MINIMUM-UNCERTAINTY GAUSSIAN WAVE PACKETS 

The standard initial minimum- uncertainty Gaussian wave packet, which gives the familiar 
time-dependent spread in Eqn. can be written in generality as 



UP) = <P(G)(p, 0) = A /4= e-" 2 ^) 2 / 2 e -W» (i 6 ) 
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where xo,Po are used to characterize the arbitrary initial central position and momentum 
values respectively. This form gives 

(P)t=Po, (P 2 )t=pl + 7T^1 and ^Pt = ^P0 = -^J= (17) 

which are, of course, consistent with the general results in Eqns. (jOJ) and (jJJ). 

The explicit form of the position-space wave function is given by Fourier transform as 

tP {G) (x,t) = — = J^= e ^-xo)/h e -« 2 (p-p ?/2 e -ip 2 mmK dj) (lg) 

V2nh V v 71 " J-oo 

which can be evaluated in closed form (using the change of variables q = p—po and standard 
integrals) to obtain 

■0(G) (Z) t) = = p ipo(x-x )/h e -ip*t/2mh p -(x-x -p t/m) 2 /2(ah) 2 (l+it/tp) 

a/ y/nahil + it/t ) 
where to = raha 2 is the spreading time. This then gives 
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(G)( 



x ,t)\ 2 = -L-e-^W^ (20) 



where 



x(t) = Xq + pot/m and /3 f = /3yi + (t/t ) 2 with (3 = ah. (21) 
This gives 

(x) t = x(t) and (x 2 ) t = [x(t)] 2 + (22) 

so that 

(Ax,) 2 = y = y( 1+ (tS) = {Axo? + {Apot/m)2 (23) 

which is the familiar textbook result, and for t = has the minimum uncertainty product 
Ax ■ Ap = h/2. 

It is easy to confirm by direct calculation that there is no initial (t = 0) x — p correlation 
(cov (x, p; 0) = 0) for this wavefunction, consistent with the lack of a term linear in t in 
Eqn. (|23|). We emphasize that such correlations do indeed develop as the wavepacket evolves 
in time, which can be seen by examining the form of either the real or imaginary parts of 
ip(G) (x, t) as shown in Fig. 1 (where we specify the model parameters used in that plot in 
the accompanying figure caption). We note that for times t > 0, the 'front end' of the wave 
packet shown there is clearly more 'wiggly' than the 'back end' (simply count the nodes 
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on either side of (x)t-) The time-dependent correlation function or covariance defined in 
Eqn. (JT3j) and correlation coefficient from Eqn. (jl4j) are easily calculated for this specific 
case to be 



cov(x,p; t) 



h ft 

2 1^ 



and 



p(x,p;t) 



t/to 



(24) 



a/i + (*Ao) 2 

which clearly expresses the increasingly positive correlation of fast (slow) momentum com- 
ponents being preferentially in the leading (trailing) ed ge o f the wave packet. We note that 



111 and 



121. 



such correlations have been discussed in Refs. 

This observation can also be described quantitatively by examining the distribution of 
kinetic energy of such a free-particle Gaussian wavepacket |l5j. In this approach, the stan- 
dard expression for the kinetic energy is rewritten using integration-by-parts in the form 



dxip*(x, t) 



d 2 ip(x, t) 
dx 2 



2m 



dx 



dip(x, t) 



dx 



which can be used to define a local kinetic energy density, T(x,t), via 



T(x,t) EE 



2m 



dtp{x, t) 



dx 



where 



(f)t 



T(x,t) dx ee T(t) 



(25) 



(26) 



As this notion is useful in systems other than for free particle states, we allow for the 
possibility that the total kinetic energy varies with time. Since this local density is clearly 
real and positive-definite, we can use it to visualize the distribution of kinetic energy (or 
'wiggliness') in any time-dependent wavefunction. We can then define similar quantities for 
the kinetic energy in the 'front' and/or 'back' halves of the wave packet, using (x)t as the 
measuring point, via 



T (+) (t) ee / T{x,t)dx 

'(*>« 



and 



T^\t) 



T(x, t) dx . 



(27) 



For the standard Gaussian wave packet in Eqn. JTHJ), the local kinetic energy density is 
given by 



T(G){x,t) = — (pI + 



2m 



2[x- x(t)]p 
a 2 h 



t/t 



[x - x(t)] 2 



+ 



l + t 2 /t 2 )\ (a 2 h) 2 (l + t 2 /t 2 



\lp(G)(x,t)\ 2 . 



The expectation value of the kinetic energy is correctly given by 



f +OQ If 1 

T ^ it) = L T ^ t)dX= 2m-{ pl+ 2^ 



(28) 



(29) 



and receives non-zero contributions from only the first and last terms in brackets in Eqn. 
since the term linear in [x — x(t)] vanishes (when integrated over all space) for symmetry 
reasons. The individual values of T^) (t) can also be calculated and are given by 



which are individually positive definite. The time-dependent fractions of the total kinetic 
energy contained in the (+)/(—) (right/left) halves of this standard wave packet are given 
by 

* (^M^) -A . (3D 



< G " ' - r<+>(t) + r ( (->(t) 2 V(2(p„q) 2 + dJ ^ttwm ' 

For the model parameters used in Fig. 1, for t = 2to this corresponds to R^/R^ = 
56%/44%, consistent with the small, but obvious, difference in the kinetic energy distribution 
seen by 'node counting'. 

Finally, this growing correlation can be exhibited in yet another way, namely through the 
Wigner quasi-probability distribution, defined by 

P w (x,p;t) = — / r(x + y,t)ij(x-y,t)e +2 ^ h dy (32) 

i r+co 

= -r / ^(p + q,t) ( p(p-q,t)e- 2 ^ h dq. (33) 

This distribution, first discussed by Wigner jl^, and reviewed extensively in the research 
~y\ and pedagogical [3| literature (and even in the context of wave packet spreading |3|)) 
is as close as one can come to a quantum phase-space distribution, and while not directly 
measurable, can still be profitably used to illustrate any x — p correlations. For the standard 
minimum-uncertainty Gaussian wavefunctions defined by Eqns. f!16|) or (|19p. one finds that 

Q 

p w ( Xj p-t) = ^- e- {p - po)2o? e -(—o-pt/mf/^ = Pw ( x _ pt / m ^ p . o) . (34) 

Contour plots of Pw{x,p;t) corresponding to the time-dependent standard Gaussian wave 
packet for two different times (t = and t = 2t ) are also shown at the bottom of Fig. 1, 
where the the elliptical contours with principal axes parallel to the x, p axes for the t — case 
are indicative of the vanishing initial correlation, while the slanted contours at later times 
are consistent with the correlations developing as described by Eqn. (|2*l*jl . (We note that 
Bohm uses a similar illustration, but discusses it only in the context of classical phase 
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space theory and Liouville's theorem.) The visualization tools used in Fig. 1 (explicit plots 
of Re[i/)(x,t)], and the Wigner function) and the distribution of kinetic energy as encoded 
in Eqns. (|50|) or ()31j). can then directly be used to examine the correlated wave packets we 
discuss in the next section. 

As a final reminder about the quantum mechanical "engineering" of model one- 
dimensional wavepackets, we recall that since an initial (po(p) is related to the time- dependent 
ip(x,t) for free-particle solutions via 

1 



ip(x,t) 



then the simple modification 



0o (p) e 



-ip I t/2mh 



Jpx/h 



dp 



(35) 



Mp)=Mp)e- ipa/h e ip2T/2mh 



(36) 



leads to the related position-space wavefunction satisfying 
1 



ip(x,t) 



Mp) e 



-ipa/h e ip T 



/2mftj e ~ip 2 t/2mh 



ipx/h 



dp = ip(x — a,t — r) (37) 



so that simple shifts in coordinate and time labels are possible, and squeezed states often 
make use of similar connections. 



IV. CORRELATED GAUSSIAN WAVE PACKETS 
A. Squeezed states 

One of the simplest modifications of a standard minimum-uncertainty Gaussian initial 
state which induces non-trivial initial correlations between position and momentum is given 
by 

(S )(p, 0) = ^= e -* 2 (P- Po f(wc)/2 e -wa_ (38) 

(A similar version of a squeezed state, but with ift(x, 0) modified, has been discussed in 
Ref. [1J|.) Because the additional C term is a simple phase, the modulus of 4>(p,t) is 
unchanged so that the expectation values of momentum, (p) Q and (p 2 )o, and the momentum- 
spread, are still given by Eqn. (|T7|) as for the standard Gaussian example. However, there 
is now an obvious coupling between the usual 'smooth' exp(— a 2 (p — po) 2 /2) term which 
describes the peak momentum values and the 'oscillatory' exp(— ipx /h) terms which dictates 
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the spatial location and spread, governed by the presence of the new C term, which leads 
to a non-zero initial x — p correlation. 

The time-dependent position-space wavefunction is obtained via Fourier transform with 
literally no more work than for the standard Gaussian and one finds 

ll} lq J x t) = \ p ip Q {x-x )/h e -ip*t/2mfh p -(x-x -p t/m) 2 /2f3 2 (l+i[C+t/t ]) /gg\ 

giving 

\^ (s) (x,t)\ 2 = -jl— e -M*>IW0 whe re b(t) = ^\ + {C + t/t,f. (40) 

Thus, the initial state in Eqn. (fHKJ) gives the same time- dependent Gaussian behavior as 
the standard case, still peaked at x = x(t), but with a spatial width shifted in time from 
t — > t + Cto- This can be understood from the results in Eqns. (JHfij) and (J37j) where the new 
C-dependent terms in Eqn. (j3*5]l give rise to effective a and r shifts given by 

a = —Ca 2 hp and r = —Ca 2 mh = —Ct . (41) 

The t shift then affects the time-dependent width, b(t), but the combined a, r shifts undo 
each other in the argument of the Gaussian exponential because they are highly correlated 
due to the form in Eqn. (JHBJ)- 

The time- dependent position expectation values are then 

(x) t = x(t) and (x 2 ) t = [x(t)] 2 + ^! (42) 



so that 



b Z.i Q 

\2 



(A Po ) 2 t 2 

,2 



(AxoY +At+ v ™' (43) 



which has a non- vanishing linear term if C ^ 0. The initial width of this packet is larger than 



for the minimal uncertainty solution by a factor of y/1 + C 2 , but has the same quadratic 
time-dependence since Ap is the same. 

One can confirm by direct calculation that 0(s)(p, 0) and ip(S)i x i 0) both do have an initial 
non- vanishing correlation leading to this form and this is also clear from plots of the initial 
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wave packet as shown in Fig. 2. We plot there an example with the same model parameters 
as in Fig. 1, but with C = —2 which leads to an anti-correlation (since C < 0) with higher 
momentum components (more wiggles) in the 'back edge' of the initial packet. This gives 
an intuitive expectation for a wave packet which initially shrinks in time, consistent with 
the result in Eqn. (|43|) . and with the plot shown in Fig. 2 for t = 2t - The parameters 
were chosen such that for this time the initial correlation has become 'undone', leading 
to something like the standard Gaussian initial state, from which point it spreads in a 
manner which is more familiar. The initial correlation is achieved, however, at the cost of 



increasing the initial uncertainty principle product by a factor of vl-j- C 2 . The complete 
time-dependent correlation coefficient from Eqn. (J 14)) is 

p(x,p;t) = ^£= +t/to) (44) 
y/l + (C + t/t ) 2 

corresponding in this case to a roughly 90% initial correlation. The required initial corre- 
lation is also clearly evident from the Wigner quasi-probability distribution for this case, 
where we find 

P w (x, p; t) = e -(P-Po) 2 ° 2 e -(x-x - P t/m-c{p- PQ )t / m f/^ _ ( 45 ) 

In this case, the initial correlation for C < shown in Fig. 2 is consistent with the desired 
anti-correlation, since the slope of the elliptical contours is negative. 

In a very similar manner, the expressions for the kinetic energy density distribution from 
Eqn. (J3*T|) are simply shifted to 

? (±) m _ T lt)(f) _ 1 ^ ( 2 \( ( Po a) \ (C + t/t ) 



R{S){t) ~ T$(t)+T$(t) 2 ± [V^J {(2( Po a) 2 + I) J yiTWTW (46) 

so that for C < 0, there is an initial asymmetry in the front/back kinetic energy distribution, 
with more 'wiggles' in the trailing half of the packet. For the C = —2 case in Fig. 2, the 
initial [t = 0) front/back asymmetry is R^/R^ = 44%/56%. 

We note that while a number of quantities (time-dependent spread in position, correlation 
coefficient, kinetic energy distribution) are simply obtained by the t — ► t + Ct shift, other 
important metrics, such as the autocorrelation function j^jj, A(t), retain basically the same 
form. 

One can imagine generating initial Gaussian states with non-zero correlations of the type 
in Eqn. (J3SJ), motivated by results obtained by the use of modern atom trapping techniques, 
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such as in Ref. [22] . In a number of such experiments, harmonically bound ions are cooled 
to essentially their ground state, after which changes in the external binding potential can 
generate various nonclassical motional states such as coherent states (by sudden shifts in the 

n 

central location of the binding potential [23]) and squeezed states (by changing the strength 
of the harmonic binding force, i.e., the spring constant). The subsequent time-development 
of Gaussian packets in such states can then lead to the desired correlated states, at which 
point the external binding potential can be suddenly removed, with free-particle propagation 
thereafter. 

As an example, the initial state in a harmonic oscillator potential of the form V(x) = 
mu 2 x 2 /2 given by 



evolves in time as ll 

ifj(x, t) = exp 

where 



imujx 2 cos (cut) 



1 



2ftsin(u;t) J y/A(t)y/n 

h 



muij3 



sin (cut) 



A(t) = f3cos(ojt) + i 

The time-dependent expectation values are then 

\A(t)\ 



exp 



and 



irruvP (x — x s (t))' 
2hsm(cut) A{t) 

Po sin (cut) 



x s (t) 



mw 



(x) t = x s (t) 



Ax t 



V2 



and 



(p)t = Po cos(ujt) 



(47) 



(4* 



(49) 



(50) 



and it is then easy to show that the time-dependent correlation of this state is given by 



cov(x,p; t) 



mio sin(a;t) cos(a;t) 



h 



mw(3 



-P 2 



(51) 



For the special case of coherent states, where (3 = Wh/mu, the correlations vanish identically 
for all times (as does the asymmetry in kinetic energy [15]), while for more general solutions, 
removing the potential at times other than integral multiples of r/2 (where r is the classical 
period) would yield an initially correlated Gaussian. 



B. Linear combinations of Gaussian solutions 



One of the simplest examples of correlated position-momentum behavior of a system, 
leading to an initial shrinking of a spatial width, can be classically modelled by two ID 
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non-interacting particles, with the faster particle placed initially behind the slower one. 
A quantum mechanical solution of the free-particle Schrodinger equation involving simple 
Gaussian forms which mimics this quasi-classical behavior, and for which all expectation val- 
ues and correlations can be evaluated in simple closed form, consists of a linear combination 
of two minimal-uncertainty Gaussian solutions of the form 



^ 2 (x, t) = N cos(0)^](:r, t) + sm(6)^(x, t) 



(52) 



where A, B correspond to two different sets of initial position and momentum parameters, 
namely (x A ,p A ) an d (xb,Pb), describes the relative weight of each component, and N is 
an overall normalization; we assume for simplicity that each component Gaussian has the 
same initial width, (3. Since each ip(G)i x ^t) is separately normalized, the value of N can be 
easily evaluated using standard Gaussian integrals with the result that 



N~ 2 = 1 + sin(2#) e 



-(x A -x B ) 2 /4l3 2 -( PA - PB ) 2 p 2 /4h 2 



cos[(x B - x A )(p B +p A )/2h] (53) 



so that if the two initial Gaussians are far apart in phase space, namely if 

(x A - x B ) 2 , (pa - Pb) 2 (3 2 



+ 



>> 1 



(54) 



A(3 2 Ah 2 

the normalization factor iV can be effectively set to unity, and all cross-terms in the evalu- 
ation of expectation values can also be neglected. 

In this limit, the various initial expectation values required for the evaluation of the 
time-dependent spread in Eqn. (J2J) are given by 

(x)q = cos 2 (8)x a + sin 2 (9)xb (55) 
(x 2 ) = co S 2 (6) (x 2 A + y) + sin2 W [ x b + y) - [™s 2 (e)x A + sm 2 (6)x B } 2 (56) 



so that 



(Ax ) 2 = [sm(26)f 



\ 2 a 
x A -x B \ B 



+ 



with a similar result for the momentum-spread, namely 

(Ap ) 2 = [sin(20)] 
The necessary initial correlation is given by 



2 (Pa~Pb\ ft 



26 2 



(xp + px) - 2(i) (p)o = 2[sin(2fl)] 2 



[x A - x B )(p A -p B ) 
4 



(57) 



(58) 



(59) 
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so that the time-dependent spread in position is given by 



(Ax t ) 2 = [sm(29)f 



%A~ %b\ ( Pa- Pb\ t 



m 



2 8 2 H 2 t 2 

4+2^' < eo ' 



In the limit we're considering, namely when \xa — %b\ » ft and/or \pa — Pb\ » h//3, the 
time-dependent width can be dominated by the quasi-classical value dictated by two well- 
separated 'lumps' of probability, and if (xa — %b) and (j>a — Pb) have opposite signs, then 
this large position spread can initially decrease in time because of the initial correlations. 
This example, while not as 'quantum mechanical' as that in Sec. II V Al does clearly and 
simply exhibit the position-momentum correlations necessary for the presence of the A term 
in Eqn. ()43|). with the 'fast one in the back, and the slow one in the front'. 

One can imagine producing linear combinations of isolated, but highly correlated, Gaus- 
sian wave packets at very different points in phase space, by invoking the dynamical time- 
evolution of bound state wave packets which leads to the phenomenon of wave packet re- 
vivals, especially fractional revivals 2^|. For the idealized case of the infinite square well 



potential 



261 ] . at t = T rev /4 (where T rev is the full revival time), an initially localized wave 



packet is 'split' into two smaller copies of the original packet, located at opposite ends of 
phase space |22[, of the form in Eqn. (|52jl . If, in this model system, the infinite wall bound- 
aries are suddenly removed at such a point in time, we then have the case considered in this 
section. 



V. CONCLUSION AND DISCUSSION 

The study of the time-dependence of the spatial width of wave packets in model sys- 
tems can produce many interesting results, a number of which are quasi-classical in origin, 
while some are explicitly quantum mechanical. Time-dependent wave packet solutions of 
the Schrodinger equation for the harmonic oscillator are easily shown to exhibit intricate 
correlated expansion/contraction of widths in position- and momentum-space j^jj and mod- 



ern experiments |22j, [2^J can probe a wide variety of such states. Even the behavior of 
otherwise free Gaussian wavepackets interacting with (or 'bouncing from') an infinite wall 

n n n 

[28], [29], j30( can lead to wave packets which temporarily shrink in size. 

While the fact that free-particle wavepackets can also exhibit initial shrinking of their 
spatial width is well-known in the physics pedagogical literature, it is perhaps not appre- 
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ciated enough in the context of introductory quantum mechanics courses because of the 
seeming lack of simple, mathematically tractable, and intuitively visualizable examples, and 
we have provided two such simple cases here. We have also emphasized the usefulness of 
several tools for the detailed analysis of the structure of quantum states as they evolve, 
namely the direct visualization of the real/imaginary part of the spatial wavefunction, the 
time-dependent spatial distribution of the kinetic energy (how the 'wiggliness' changes in 
time), and the Wigner quasi-probability distribution all of which provide insight into the 
correlated x — p structure of quantum states. 



APPENDIX A: TIME-DEPENDENT Ax t FOR THE CASE OF UNIFORM AC- 
CELERATION 



Many of the results discussed here for the time-dependent widths of localized quantum 
wavepackets for the free-particle case can be carried over to the situation of a particle under- 
going uniform acceleration, governed by the Hamiltonian H = p 2 /2m — Fx, corresponding 
to a constant force, +F, to the right. General expressions for the time-dependent values of 
expectation values of powers of both x and p can be obtained using Eqn. © (following the 
method of Styer [8( for example) and one obtains the results 

(p)t = Ft+(p) and (p 2 ) t = FH 2 + 2F(p) t + (p 2 ) (Al) 

which imply that (Ap f ) 2 = (Ap ) 2 . For position we have the corresponding results 

Ft 2 (p)nt , . 

*>t = 7T" + — + *>o (A2) 
2m m 

a FH 4 F(p) t 3 (p 2 ) t 2 F(x) t 2 (xp + px) t 2 

\ x f = t~^. + — t~ + — r~ + + + x A3 

which actually combine to give the same expression for (Axt) 2 as in Eqn. (J2J). Thus, the 
model systems we have discussed here can also be used as examples of correlated wave 
packets in the related accelerating particle case. 

The expressions above can also be derived using an approach similar to that followed in 
Sec. |nj namely by using the most general form for the time-dependent momentum-space 
wavefunction, <j>(p,t). In that case, the time-dependent Schrodinger equation in p-space is 
written in the form 



op 



<f>(p,t)=ih^<f>(p,t) (A4) 
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which has a general solution 

0(p, t) =Mp~ Ft) S p - Ft)3 - p3]/QmFh (A5) 

where </>o(p) = </>(p, 0) is still the initial wavefunction. The t-dependent expectation values 
in Eqn. (jAlj) . ()A2|) and ()A3|) can then be obtained as in Sec. |H]in terms of the t — results, 
just as for Eqns. (jgj), (171). ©, and (fTHIl. 
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Figure Captions 



Fig. 1. The time-development of a standard Gaussian wavepacket, tf)^(x,t), described by 
Eqn. with the modulus (solid) and real part (dashed) shown for t = and 

t = 2t on the top plot. Contour plots of the Wigner function, from Eqn. (}3"4*j) . for the 
same two times are shown at the bottom, corresponding to contours which are 70%, 
30%, and 10% of the peak or central value. The model parameters used in this plot 
are H = m = 1 and (3 = 2, which give to = m/3 2 /h = 4, along with p = 4 and Xq = —4 
for the initial packet. 

Fig. 2. Same as Fig. 1, but for the correlated squeezed state in Eqn. (J3~%j) or (j3Uj) . with the 
model parameters used in Fig. 1. For this case, we use C = — 2, so that for t = 2t 
the initial correlations are 'undone'. 
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contour plot of P w (x,p;t) vs. (x,p) 
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